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In previous works Suppes and de Barros used a pure particle model to 
derive interference effects, where individual photons have well-defined tra- 
jectories, and hence no wave properties. In the present paper we extend 
that description to account for the Casimir effect. We consider that the 
linear momentum ^ ^Tik of the vacuum state in quantum electrodynamics 
corresponds to the linear momentum of virtual photons. The Casimir effect, 
in the cases of two parallel plates and the solid ball, is explained in terms of 
the pressure caused by the photons. Contrary to quantum electrodynamics, 
we assume a finite number of virtual photons. 
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1 Introduction 

Suppes and de Barros (1994a, 1994b, 1996) began a foundational analysis on 
diffraction of light which formulated a probabilistic theory of photons with 
well-defined photon trajectories and without wave properties. The wave 
properties come from the expectation density of the photons. The photons 
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are also regarded as virtual, because they are not directly observable, in- 
cluding their anihilation of each other (see assumption (5) below). What 
can be detected is the interaction with matter. The meaning of virtual used 
here is not the same as in quantum electrodynamics (QED). In summary, 
our assumptions are: 



• Photons are emitted by harmonically oscillating sources; 

• They have definite trajectories; 

• They have a probability of being scattered at a slit; 

• Absorbers, like sources, are periodic; 

• Photons have positive and negative states (+-photons and — photons) 
which locally interfere, when being absorbed; 

• Photons change their states when reflected by a perfect conductor. 



The expected density of ib-photons emitted at t in the interval dt is given 

by 

s±{t) = cos iVt), (1) 



where lo is the frequency of a harmonically oscillating source. As is a constant 
determined by the source, and t is time. If a photon is emitted at t', < 
t' < t, then at time t the photon has traveled (with speed c) a distance r, 
where 

t-t'='^. (2) 

The conditional space-time expectation density of ib-photons for a spheri- 
cally symmetric source with given periodicity uj is: 

h±{t,r\u;) = -^(^l±cosij(^t-'^^y (3) 

where A is a real constant. 

The scalar field defined in terms of the expectation density h±{t, r\Lo) is 

h+ - h- 
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where £o is a scalar physical constant. Using (|3|), @ may be rewritten for 
a spherically symmetric source as: 



^ = ^o\/^cos..(t-0. (5) 

In the present paper we construct a corpuscular model for the Casimir 
effect, following the ideas ever proposed by Suppes and de Barros (1994a, 
1994b, 1996). First, we study the case of the parallel plates, and then the 
solid ball. 



2 Quantum Vacuum 

A pure particle theory must postulate properties of the quantum vacuum 
if standard field-theoretic methods of computing the Casimir effect, and 
similar phenomena such as the Lamb shift, are to be closely approximated. 

In QED the vacuum state has an energy ^hoj and a linear momentum 
^Kk. We consider, as a first postulate, that ^Tik corresponds to the linear 
momentum of one virtual photon. Our second postulate is that we have 
a distribution function of k, f{k), which is a probability density for the 
distribution of photons with respect to k. Our third postulate establishes 
that both +-photons and — photons contribute to the delivery of linear 
momentum on a reflective surface. So, the conditional expectation density 
of photons, given k, that strike a point on the conductor surface is 

h{t,rs\k) = h+{t,rs\k)+h^{t,rs\k), (6) 

where rs is a surface point. 



3 Parallel Plates 

We consider here the case of two perfectly conducting parallel plates, stand- 
ing face to face in vacuum at a distance d much smaller than their lateral 
extensions. It is well known that such plates attract each other with a force 
per unit area (pressure) due to the vacuum energy, as predicted by Casimir 
(1948), given by 

P = -^ (7) 
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Usually, such an attraction is explained in terms of the vacuum field. We use 
a random distribution of oscillating sources of photons, in the vacuum, which 
do not interfere with each other, to derive in our conceptual framework. 

The photons outside the plates that strike such surfaces act to push the 
plates together, while reflections of the photons confined between the plates 
push them apart. This idea is proposed by Milonni, Cook, and Goggin 
(1988) and also presented in (Milonni, 1994), but not actually developed 
from a pure particle viewpoint. 

The photons that we are considering must satisfy a probability den- 
sity f{k) > 0. Rather than assume an explicit expression for f(k) (which 
requires some assumptions about the virtual photons), we prefer to state 
some properties that f{k) must satisfy: 

(i) f{^)dkxdkydkz = 1, and the mean and the variance of f{k) 
are finite; 

(ii) There exists a constant H such that h{t,rs\k) < H. 

(iii) h(t,rs\k)f{k)\k=Q = 1, and all derivatives of this expression vanish at 

fc, = 0. 

From (i)~(iii) and assumptions made earlier, we may infer, contrary to 
a standard result of QED, that in our theory the number of virtual pho- 
tons is finite for any bounded region of space-time. We also infer that 
h{t,rs\k)f{k)\k=oo = Oi which is intuitively an expected property of a cut- 
off function. We note that h{t,rs\k) and f{k) are not identified as specific 
functions. We have generalized from standard cutoff functions, such as an 
exponential function, to give reasonable sufficient conditions that many dif- 
ferent functions satisfy. We do not know enough about the quantum vacuum 
to derive a particular choice. 

We divide the xyz space into parallelepipeds of sides L^-, Ly, and L^, as 
in the usual description of QED. So, all k^, ky, and kz must assume discrete 
values, as is explained in the next paragraphs. 

We note that when reflected a photon changes its state from positive to 
negative and vice versa. This single change for perfect conductors implies 
that the defined scalar field, given by (Q), vanishes at the refiecting surface. 
For further details see (Suppes and de Barros, 1996). So, according to (^) 
and recalling that k = uj/c, we have at the wall: 

cos (ivt — Lo- ) = cos {ujt — kr) = 0. (8) 
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If we set ujt = tt/2, which corresponds to a convenient choice for the origin 
of time, it is easy to see that the values of k^, ky, and kz that satisfy the 
boundary condition in x = L^, y = Ly, and z = Lz are: 

^^'i^'^ = ^ (9) 

But the condition that the scalar electric field vanishes at the surface of 
the reflecting walls is not sufficient to explain the periodicity given by (^). A 
natural question arises: what about the photons with linear momenta that 
do not satisfy @)1 We recall that reflectors, like absorbers (Suppes and 
de Barros, 1994b), behave periodically, since the photons are continuously 
hitting the plates. Thus, the probability of reflecting a photon is given by: 

p = C(l + cos(a;t + V)), (10) 

where ijj is a. certain phase. If p = 0, then there is no interaction with the 
plates, which means that no momentum is delivered to it. 

As an example, consider the first strike of a photon on a plate perpen- 
dicular to the z axis. Such a surface is not oscillating before the strike. But 
after reflection, the wall oscillates with the same frequency uj associated to 
the linear momentum k = ujjc. The particle reflects on the other wall and 
returns to the first wall with a phase 2Lzkz- But we must have 2Lzkz = 2mr, 
from (|l^), if the particle is to be reflected again on its return to the first 
wall. Obviously, cos(a;t — 2Lzkz) = cos{ujt) if and only if 2Lzkz = 2mr. 

If we consider L^^y^z very large compared with any physical dimensions 
of interest, we can assume that the kx,y,z approach a continuum. This is 
what holds for photons outside the plates. 

Now we start to derive the pressure on the plates. We begin with the 
inward pressure. The expected number of photons that strike the area dS 
of one of the plates, within the time interval dt is 

h{t,rs\k)f{k)^dkxdkydkz cosjcdtdS, (11) 

where 7 is the angle of incidence of the photons on the plate with respect 
to the normal of the surface, i.e., C0S7 = kz/k, where k = ^k]. + k'^ + k'l. 
Thus, the element of volume that we are taking into account is cos jcdtdS. 
The factor ^ is justified by (^), since outside the plates we approach the 
continuum as a limit. 



5 



The momentum dehvered to the plate by a single reflected photon is 
equal to the negative of the change in the momentum of the photon. In 
other words the momentum is equal to 2^fikz, if we consider the plate per- 
pendicular to the z component of the xyz system of coordinates. Therefore, 
the expected linear momentum transfered to an area dS on the plate during 
the time interval dt is 

-^-^h{t,rs\k)f{k)dk^dkydkzcdtdS. (12) 



The force on the plate is obtained by dividing (|T2|) by dt. The pressure 
is obtained by dividing the force by dS. We denote the inward pressure as 
Pin and the outward pressure as Pout- Hence: 

dPin — q , -dkxdkydkz- (13) 

TT'^ /Zc-2 ^ i.2 I 1,2 

Integrating over momentum: 

p.„ = 5^ r r r ,,jii<Mm. 

^^^0 Jo 'Jo Jkl + kl + kl 



The equation given above is identical to a result due to Milonni, Cook and 
Goggin (1988), if we consider that h{t,rs\k)f{k) has the role of the usual 
cutoff function. 

To obtain the expression of the outward pressure we use similar argu- 
ments. But now, because of the small distance d between the plates, we must 
take into account the periodicity given in (P), at least for the z component. 
Hence' 

p r r Ht,rs\k)f{k){^f 

Pout = dkx dky — . (15) 

-^'d^Jo Jo .^kl + kl + (^)2 

We note that it follows from (i)~(iii) that Pin and Pout are both finite. 
The resultant pressure is given by: 

Pout Pin — 



■k'^Tic 2 /"°° T h{t,rs\x^u)f{^/x'■ 



— ^ (-OO 

■ I dx 



Jo 



n 



2 



Tr'^hc f°° 2 f°° 1 h{t,rs\x,u)f{Vx + u 



2\ 



Jo io ' Vx + .^ ' 
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where, by change of variables, f{k) = f{\/x + v?), x = x'"^ = ^ — h 
u = kz^, 9 = tan^^^, and dkxdky = x'dx'd9^. The expres- 
sion h{t, rs\x, u)f{Vx + n^) corresponds to a cutoff function. In our model 
h{t,rs\k) is bounded and f{Vx + u^) has the physical interpretation of a 
probability density of the frequencies of the photons. 

Frequently it is assumed that the cutoff function has the property of 
going to zero as k approaches infinity and going to one when k approaches 
zero. This is justified physically with the hypothesis that the conductivity 
of the reflecting conductors decreases to zero as the frequency gets high. 
Since h{t,rs\k) is bounded, it is easy to see that the product h{t,rs\k)f{k) 
must assume a similar role with respect to the cutoff, from a mathematical 
standpoint. According to the assumptions that we made about h(t, rs\k) and 
f{k), all the properties of a cutoff function are satisfied for h{t,rs\k)f{k). 

If we consider: 



N 2 f°° ^ h{t,rs\x,u)f{Vx + u^) 

F{u) = u / dx , , , (17) 

JO yx + w^ 

it is clear that the Euler-MacLaurin summation formula (Abramowitz and 
Stegun, 1971) may be applied to (|T^). Contrary to the standard QED 
treatment, we apply this formula to a convergent rather than a divergent 
series. Therefore, the factor that is multiplying in (IC) may be written 



as: 

°° /-oo 1 1 1 

E i^H - X duF{u) = --F(O) - -F'(0) + 7^i^'"(0)... (18) 

for lim^i^oo = 0, since Yl,^=i is finite and so lim^^oo F{n) = 0. 
We note that F(0) = 0, F'(0) = 0, F'"(0) = -12/i(t, rs|0), and all higher 
derivatives F(")(0), where n is odd, vanish in accordance with assumption 
(iii), as is shown in the appendix. Since by (iii) h{t,rs\Q)f{^) = 1: 

P _p -_^!^ ^^q^ 

Equation ( p^ ) is identical to (^), which completes our derivation of the 
Casimir effect for parallel plates. 

4 The Solid Ball 

Many authors have considered the case of the Casimir force for solid balls 
and cavities like spheres, hemispheres, and spheroids. Balian and Duplantier 
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(1977) developed a method that estabhshes an expansion for the Green 
functions describing electromagnetic waves in the presence of a perfectly 
conducting boundary. Later, they applied their method to the study of 
the Casimir free energy of the electromagnetic field in regions bounded by 
thin perfect conductors with arbitrary shape (Balian and Duplantier, 1978). 
Brevik and Einevoll (1988) used Schwinger's source theory to establish the 
Casimir surface force in the case of a solid ball, considering e{uj)fi{Lo) = 1, 
where e{uj) is the spectral permittivity and /u(w) is the spectral permeabil- 
ity. In 1990, Brevik and Sollie (1990) calculated the Casimir surface force 
on a spherical shell, assuming the same condition e{uj) = 1. Barton 
(1991a, 1991b) uses standard statistical and quantum physics to analyze 
the fluctuations of the Casimir stress exerted on a flat perfect conductor by 
the vacuum electromagnetic fields in adjacent space. Eberlein (1992) makes 
an extension of Barton's work, calculating the mean-square forces acting on 
spheres and hemispheres of variable sizes. 

We adopt here the same corpuscular model presented in last section, to 
the case of a solid ball of radius a surrounded by vacuum. Consider a ball, 
centered at the origin of a spherical coordinate system {p,(p,9). At each 
point of the surface of the ball, we define a Cartesian system of coordinates, 
with axes _L, ||i, and ||2. -L is the normal to the surface, and ||i and ||2 are 
tangent to the sphere. 

As in the case of the parallel plates, we assume a distribution function 
f{k) satisfying the same properties assumed in the last section. The ex- 
pected number of photons that strike the area dS = o? sin ^p dip dO on the 
surface of the solid ball, within the time interval dt is 



where 7 is an angle of incidence of the photons with respect to the normal 



As in the case of the plates, the momentum delivered to the ball by a 



h{t, rs\k)f{k)dk_\_dk\\^dk\\^cdt 00370^ sm(pd(pd6. 



(20) 




single reflecting photon is 



2-^hk±. The linear momentum on the ball is 



hk'^. 



h{t,rs\k)f{k) 

TT^k 



dk^dk^^-^dk^^^cdt a^sm(pdipd9. 



(21) 



The force is obtained by dividing the expression above by dt: 
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Integrating: 



it' Ja Ja "'A '\/fc2 +B +B 



2 



By the arguments used earher, the force F is finite. 

Our result depends explicitly on h(t,rs\k)f{k) in (p3|), which has, as in 
the case of the parallel plates, a role similar to a cutoff. This is a consequence 
of the geometry of the problem. Brevik and Einevoll (1988) obtained another 
expression for the Casimir force in the case of a solid ball, which directly 
depends on a typical value (3 x 10^ see" ^) for the cutoff frequency. 
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Appendix 

It follows from (^) that 

F{u) = + ^ r dxg{^/^^)\{xW)-^/\ (24) 

"'0 yx + v? Jo 2 

where g{\/ x + u^) = h{t, rs\x, u)f{V x + n^). If we make y = {x + v})^!"^: 

POO 

F{u)=2u^ dyg{y). (25) 

J u 

Before the evaluation of the derivatives of F(u), we must observe that 
d 

-j^j^ dyg{y) = -g{u), (26) 

since ^(oo) = 0. Hence: 

/■oo 

F'{u)=Au dyg{y)-2u^g{u), (27) 

J u 

F"'{u) = -I2g{u) - 12ug'{u) - 2u^g"{u), (28) 
and all higher derivatives ^(^n+i) vanish at u = if the even derivatives 
of g{u) vanish at the same point, which is assumed in (iii). 
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